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Abstract —This paper proposes a novel approach for detecting 
the topology of distribution networks based on the analysis of time 
series measurements. The time-based analysis approach draws on 
data from high-precision phasor measurement units (PMUs or 
synchrophasors) for distribution systems. A key fact is that time- 
series data taken from a dynamic system show specific patterns 
regarding state transitions such as opening or closing switches, 
as a kind of signature from each topology change. The proposed 
algorithm here is based on the comparison of the actual signature 
of a recent state transition against a library of signatures derived 
from topology simulations. The IEEE 33-bus model is used for 
initial algorithm validation. 


I. Introduction 

Different tools have been developed and implemented to 
monitor distribution network behavior with more detailed and 
temporal information, such as SCADA, smart meters and line 
sensors. Creating observability out of disjointed data streams 
still remains a challenge, though. Given the present monitoring 
technologies, more, better and faster data from behind the 
substation will be needed to realize smart distribution networks 
m. The cost of monitoring systems in distribution networks 
remains a barrier to equip all nodes with measurement devices. 
To some extent, a capable Distribution System State Estimation 
can compensate for the lack of direct sensor data to support 
observability. However, switches status errors can easily be 
misinterpreted as analog measurement errors (e.g. voltage or 
current readings). Thus topology detection is an important 
enabling component for state estimation as well as a host of 
other operation and control functions based on knowledge of 
the system operating states in real-time. 

Topology detection in distribution networks is not a trivial 
problem. Firstly, switches may not reliably communicate their 
status to the distribution operator, so the topology can only 
be determined by sending crews into the field. Moreover, the 
reported switch status may be faulty due to switch malfunction 
and unreported maintenance crew manipulation. Integration of 
distributed energy resources, electric vehicles and controllable 
loads add more dynamics to the distribution networks that 
may lead to more frequent protection and switching actions. A 
survey of utilities experts done by the authors shows on average 
5 to 10 switching actions happens under an urban distribution 
substation m. Furthermore, knowledge of the correct, updated 
topology is essential for safety and service restoration after 
outages ll^. 
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Most of literature on topology detection are based on state 
estimator results and measurement matching with different 
topologies. In H authors propose a state estimation algorithm 
that incorporates switching device status as additional state 
variables. A normalized residual test is used to identify the 
best estimate of the topology. State estimation based algorithms 
are limited to state estimator accuracy. They are also sensitive 
to measurement device placement. In 13, the authors provide 
a tool for choosing sensor placement for topology detection. 
Given a particular placement of sensors, the tool reveals the 
confidence level at which the status of switching devices can 
be detected. Authors in |6l and Q are focused on estimating 
the impedance at the feeder level. However, even a perfect 
identification of network impedance cannot always guarantee 
the correct topology, since multiple topologies could present 
very similar impedances. 

The work in used power flow analysis for matching 
substation loading and aggregated household meter load data 
from network connectivity modeling. They assumed metering 
load data are time synchronized with a measurement device 
on each transformer. The assumption is still far from an actual 
load metering system. Moreover, convergence of the proposed 
optimization is sensitive to bad data. Voltage measurement 
cross correlation for house meters data is presented in i). 
However, residential meters do not usually provide voltage 
measurements for utility operations. Voltage correlation-based 
methods in balanced feeders, feeders with PV resources, 
and feeders with inaccurate GIS models can be error prone. 
Moreover, voltage measurement in meters are hourly or 15 
minute average values. Voltage average values introduce addi¬ 
tional errors in voltage regression-based methods for topology 
detection. Most of the proposed methods in literature are post¬ 
processing methods which depend on correct execution of state 
estimation or power flow. 

In this paper, a novel approach for topology detection is 
proposed based on time series analysis of phasor measurement 
unit (PMU) data. This approach is inspired by high-precision 
phasor measurement units for distribution systems micro- 
synchrophasors (/r-PMU), which the authors are involved in 
implementing nni. The main idea derives from the fact that 
time-series data from a dynamic system show specific patterns 
regarding system state transitions, a kind of signature left from 
each topology change. The algorithm is based on the compar¬ 
ison of the trend vector, built from system observations, with 
a library of signatures derived from the possible topologies 
transitions. The topology detection results are impacted by load 
uncertainty and measurement device accuracy. Therefore, the 
conduced analysis takes load dynamics and measurement error 
into account. A set of actual secondly load measurements are 
used for a number of residential customers in the United States. 
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The derived statistical load model is applied to the topology 
detection scenarios to validate the proposed algorithm. The 
topology detection accuracy is also dependent on p-PMU 
placement. We proposed a p-PMU placement approach for 
topology detection application. The analysis shows that the 
topology detection algorithm converges even with limited 
measurement devices. 


We assume that the PMUs deployed in the grid take mea¬ 
surements at the frequency /. We denote with the unique 
topology whose switches status is described by a. Its bus 
admittance matrix is defined as 


(y.hk 


if j = fc 

—Yjk, otherwise 


( 1 ) 


II. Distribution Grid Model 

This section presents the distribution network model and its 
related notations. Given a matrix W, we denote its transpose 
by and its conjugate transpose by W*. We denote the 
matrices of the real and of the imaginary part of W by 3fi(kP) 
and by 5 (VP), respectively. We denote the entry of W that 
belongs to the j-th row and to the fc-th column by Wjk- 

Given a vector v, vj will denote its j-th entry, while v-j the 
subvector of v, in which the j-th entry has been eliminated. 
We denote its complex conjugate by v. Given two vectors v 
and w, we denote by {v, w) their inner product v*w. We define 
the column vector of all ones by 1. 

We associate with the electric grid the directed graph Q = 
{V,8) and the sets S and V, where 

• V is the set of nodes (the buses), with cardinality n; 

• £■ is the set of edges (the electrical lines connecting the 
buses), with cardinality w\ 

• S is the set of switches (or breakers) deployed in the 
electrical grid, with cardinality r 

• 7^ is the set of nodes endowed with voltage phasor 
measurement units (PMUs), with cardinality p 

Let A € {0,±1}“'^" be the incidence matrix of the graph Q, 

A=[a{ ... al]'^ 

where aj is the j-th row of A. The elements of A are 
all zeroes except for the entries associated with the nodes 
connected by the j-th edge, for which they are equal to -fl 
or —1 respectively. In this study, we consider the steady state 
behavior of the system, when all voltages and currents are 
sinusoidal signals waving at the same frequency loq. Thus, 
they can be expressed via a complex number whose magnitude 
corresponds to the signal root-mean-square value, and whose 
phase corresponds to the phase of the signal with respect to an 
arbitrary global reference. Therefore, x represents the signal 
x{t) = |a;|-\/2sin(a;ofThe system state is described by 
the following quantities: 

• u G C”, where is the grid voltage at node v, 

• i € C", where iy is the current injected at node v; 

• s = p + iq € C”, where s„, and are the complex, 
the active and the reactive power injected at node v, 
respectively; 

• (T G {0,1}'’, where is the staus of the breaker v: 
(T„ = 0 if the switch v is open, (t„ = 1 if the switch v 
is closed; 

• y G CP, where yy is the grid voltage measured by the 
sensor v 

• the trend vector G C^, defined as the difference 

between phasorial voltages taken at the two time istant 
ti and t 2 . i.e. S(ti,t 2 ) = u{ti) - u{t 2 )- 


where Yjk is the admittance of the branch connecting bus j 
and bus k and where we neglect the shunt admittances. From 
Q, it can be seen that Y„ is symmetric and satisfies 

Y,1 = 0 , ( 2 ) 

i.e. 1 belongs to the Kernel of Y Furthermore, it can be 
shown that if Q, the graph associated to the electrical grid, is 
connected, then the kernel of Y^ has dimension 1. 

We model the substation as an ideal sinusoidal voltage 
source (slack bus) at the distribution network nominal voltage 
Ljv- We assume, without loss of generality, that C/at is a real 
number. We model all nodes but the substation as constant 
power devices, or P-Q buses. In the following, ui,si,pi,qi 
will denote the voltage of the substation, and its complex, its 
active and its reactive power injected, respectively. The system 
state satisfies the following equations 


i = Y cU 


(3) 

e 

II 


(4) 

llyiy - Py “h i^y 


(5) 


The following Lemma El introduces a particular, very 
useful pseudo inverse of Y 

Lemma 1: There exists a unique symmetric, positive 
semidefinite matrix G C"^" such that 


(X^Y^ = I-lel 
IX^ei =0. 


(6) 


Applying Lemma [T] from 01 and we can express the 
voltages as a function of the currents and of the nominal 
voltage Uf^: 

u = X^i + 1 Un (7) 

The following proposition (El) provides an approximation of 
the relationship between voltages and powers. 

Proposition 1: Consider the physical model described by 
the set of nonlinear equations 0, 0, 0 and 0- Node 
voltages then satisfy 

it=Ljvl+py— X„s + o (——) (8) 

Un \^nj 

(the little -0 notation means that lim( 7 jy_>.oo — 0)- 

Equation 0 is derived basically from a first order Taylor 
expansion w.r.t. the nominal voltage Un of the relation among 
powers and voltages. Even if the solutions of 0 could be 
not very precise, and then its application to the power flow 
computation could be not satisfactory, it has been already used 
with success in state estimation El, VoltWar optimization 
ifTSl . and the optimal power flow problem iflTll . 
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III. Identification of Switching Actions 

The basic idea behind our proposed approach is that changes 
of the breakers status will create specific signatures in the 
grid voltages. In order to develop the theoretical base for the 
proposed algorithm, we make the following assumptions. 

Assumption 1: We assume that all the lines have the 
same resistance over reactance ratio. Therefore, = 

Assumption 2: We assume that only one switch can change 
its status at each time. 

Assumption 3: We assume that the graph associated with 
the grid is always connected, i.e. that there are no admissible 
state in which any portion of the grid remains disconnected. 
Assumption will be relaxed in Section [VIII| in order to test 
the algorithm in a more realistic scenario. However, it allows 
us to write the bus admittance matrix in a simpler way. In fact, 
due also to Y„ symmetry, we can write 

Y,=^[Y,) + i^{Y,) = {l + ia)U^RU* (9) 

where Yir is a diagonal matrix whose diagonal entries are 

the non-zero eigenvalues of 3fi(Ycr)j and U is an orthonormal 

matrix that includes all the associated eigenvector. From (|^, 

it can be shown that U spans the image of / — 

i.e. the space orthogonal to 1. The matrix can be written 

as 

= (l + ia)-iAC/S^i[/*A^ (10) 

with A = (/ — lef^). Assumption 0 is instead reasonable 
for the proposed algorithm framework it works in the scale 
of some seconds, while typically the switches are electro¬ 
mechanical devices and their actions are not simultaneous. 
Finally, Assumption is always satisfied during the normal 
operation. We introduce the main idea that underlies our 
algorithm with the following example. 

Example 1: Assume that at time t — \ the switches status 
is described by a{t — 1) = cji, resulting in the topology 
with bus admittance matrix Y. Applying Proposition H and 
neglecting the infinitesimal term, we can express the voltages 
as 

u{t-l) = X„,^ + lUN (11) 

Un 

At time t the £-th switch, that was previously open, changes 
its status. Let the new status be described by a{t) = a-i, 
associated with the topology Since we are basically 

adding to the graph representing the grid the edge on which 
switch £ is placed, we can write 

Y,^=Y„^+YiaiaJ (12) 

where >£ is the admittance of the line, and the elements of ai 
is the £-th row of the adjacency matrix associated with T^s- 
Since is orthogonal to 1, there exists be such that Ube = ae- 
This allow us to write 


The trend vector 6{t, t — 1) can thus be written as 

= (15) 

Un 

Example 2: Consider now the opposite situation; at time 
t — 1 the f-th switch is closed and it is opened at time t, i.e. 
a{t — 1) = (72, o'{t) = ai- In this case we have that 

Sit,t-l) = iX,,-X^,)-^ (16) 

Un 

We can observe that when there is a switching action, the 
voltage profile varies in a way determined by the particular 
transition from a topology to another. Furthermore, opening 
or closing a determined breaker makes the system vary in the 
opposite way, e.g. compare ( [T5| ) with ( [T6| l. This is basically 
due to the fact that 

Given a breaker status a associated with the topology Y^ = 
{l + ia)U13iiU*, and fixed an edge £ endowed with a breaker, 
associated with the row ae in the incidence matrix, we define 
the signature matrix 

- t/(Efl + ^iYe)bebJ)-^U* (17) 

Exploiting the signature matrix we can write 

<5(f,f-l)=A$,(t)_,A^^ (18) 

Un 

if the switch £ has been closed, else 

<5(f,f-l) = -A$,(*)_,A^^ (19) 

Un 

if the switch £ has been opened. The following Proposition 
show a characteristic of the signature matrix that is fundamen¬ 
tal for the development of our topology detection algorithm. 

Proposition 2: For every transition from the state described 
by a{t — 1) to the one described by a{t), due to a action of 
the switch £, is a rank one matrix. 

Proof: Following the same reasoning of Example [T] 
and exploiting K. Miller Lemma (113), after some simple 
computations we can write 

^a(t).,=pUYf,^bebjYf,^U* (20) 


^ 1 -f Tr(3fi(r^)&,6f Efl) 

It’s trivial to see that ^w{t)_e is a rank one matrix with 
eigenvector 

( 21 ) 

associated with the non-zero eigenvalue 


Y^, = (1 -T ia)U{T.R + 'iR{Yi)bthJ)U* 

= (1 -f ia)“^A[/(Efl -T iSi{Ye)bebJ)-^U* (13) 
The voltages are 

uit) = X,^^ + lUN (14) 


=p||C/E-i6,f 

Thus, we have 




( 22 ) 
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The trend vector — 1) represents how the opening or the 
closure of a switch spreads on the voltages profile. Thanks to 
Proposition we can write it as 


S{t, t — 1) 




^ I 


N 


^9<y{t)_e 


(23) 


from which we see that 6{t, f — 1) oc i-6-> every pat¬ 

tern that appears on the voltage profile due to switching actions 
is proportional to Aga-(t)_i, irrespective of other variables such 
as voltages u and loads s that describe the network operating 
state at the time. Thus, 9a(t)_i can be seen as the particular 
signature of the switch action. This fact is the cornerstone for 
the topology detection algorithm in this paper. 

Remark 1: The opening and the closure of the switch £, 
once the other switches status is fixed, share the same signature 
9 <T{t)_t _7 thus in principle they are indistinguishable. Without 
any other information, the trend vector can just identify witch 
switch has changed its status. 


If the projection is greater than the threshold, the associated 
switch is selected and the topology change time is detected. 
If there is no switches action, the trend vector will be zero 
as all the 0 ^( 4 )and the algorithm will not reveal any 
topology transition. Thus, the projection value is used by the 
algorithm to detect the change time too, differently of what 
proposed in ifl^ . where we used the norm of a matrix built 
by measurements (the trend matrix). The new approach is more 
reliable in the realistic case. With a slight abuse of notation, 
we will say that the maximizer of C, denoted by maxC, is 
the switches status cr such that \cr\-i = [a]i = 1 if 

= 0 or vice-versa [a]^ = 0 if a{t)-e = 1 and 
its the maximum element in C. 

We tacitly assumed so far all the buses endowed with a 
PMU, just to show the main idea. But this is not a realistic 
scenario for a distribution network. Now we show how the 
former approach can be generalized in presence of limited 
information, in which we are allowed to take only a few 
voltage measures: 


IV. Topology Detection Algorithm 

If we assume the distribution network physical infrastructure 
known, i.e. conductor impedances and switch locations, we 
can build a library C in which we collect all the normalized 
products between A and the eigenvectors ( |2T] i for all possible 
breaker actions 


^ = {9w_ ^ : w satisfies Assumption (24) 


where 

Apcr(t)f 


(25) 


It is natural to compare at each time the trend vector 
with the elements in the library to identify if and which switch 
changed its status. As stated in Remark [T] if we want to 
identify which is the current topology we need additional 
information, i.e. the knowledge of the topology before the 
transition. In that case, we could compare the trend vector 
with a restricted portion of the library £, since there are only 
r possible transitions, each of one caused by the action of one 
of the r breakers. As a consequence, we can compare S{t, t—1) 
with the particular library 


^a(i-i) = {9<T{t)_i ■ cr{t)-i = a{t - l)-^} (26) 


that is peculiar of the state a{t — 1) before the transi¬ 
tion. The comparison is made by projecting the normalized 
measurements-based trend vector ||^|*’*~]^||| onto the topology 
library The projection is performed with the inner 

product, and it allows us to obtain for each vector in £o.(t-i) 





(27) 


If Ccr(t)_f — it means that 5 is spanned by 9a(t)_i and 
then that the switch i changed its status. Because of the 
approximation I®., the projection will never be exactly one. 
Therefore, we will use a heuristic threshold, called min_proj, 
based on our numerous simulations to select the right breaker. 


y = Ivu 

= /pX,(t)^ + lC/jv (28) 

Un 

where I-p G [0, 1]p^" is a matrix that select the entries of u 
where a PMU is placed. In that case, the trend vector becomes 

S{ti,t 2 ) = yih) - y{t 2 ) (29) 


The library vectors and their dimension change too. In fact 
one can easily show, using ( [28] l and retracing •HI and ( |20| i 
that (|25|l becomes 


_ Iv^ga{t)_e 

" UvAg^wJ 


(30) 


The topology detection algorithm with limited measurements 
is stated in Algorithm [T] 


Algorithm 1 Topology Changes Detection 

Require: At each time t, a{t — 1), min_proj = 0.98 
I: a{t) G- a(t — 1) 

2: each PMU at each node j record voltage phasor measure¬ 
ments yj{t) 

3: the algorithm builds the trend vector 6{t,t — 1) 

4: the algorithm projects S{t,t — 1) in the library C-p^^rit) 
obtaining the set of values 



^y9a(t)_, 


1 9cr{t)-i 



5: if maxC > min_proj then 
6: a(t) G- argmaxC 

7: end if 


V. PMUS PLACEMENT 

When we have a limited number of sensors to be deployed 
in the distribution grid, the first requirement to be satisfied 
is the system observability, i.e. the algorithm capability of 
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detecting every topology transition. Since our algorithm is 
basically based on the comparison between trend vectors and 
the library vectors, the trivial condition for the observability of 
the network is that each vector of the library is not proportional 
to any of the others. The former property is equivalent to the 
following condition. 

Lemma 2: Given the set of nodes endowed with PMUs V, 
let, with a slight abuse of notation, the juxtaposition of the 
library vectors be denoted by C. Let C-p = IpC. Then if 

{\C*pCv\)uv<iyu^v (31) 

the switch that changes its status can be identified. 

Notice that the element {CpCp)uv is simply the projection of 
9[a{t)]-^ onto g[cr(t)]_„. If they are note purely proportional, 
trivially the magnitude of their inner product is smaller than 
one. Condition iB can be too restrictive, if we know the 
switches status before the topology change. In that case, in 
fact, we can just check if each vector of the particular library 
is not proportional to any of the others. The former property 
is equivalent to the following condition. 

Lemma 3: Let the switches status be cr, and let it be known. 
Given the set of nodes endowed with PMUs V, let, with a slight 
abuse of notation, the juxtaposition of the particulary library 
vectors be denoted by Ca- Let Cp^a = IpCa- Then if 

< l,Vu^t;,V(T (32) 

each switch action can be identihed. 

Lemmas and can be used to infer, given the electrical 
grid model, which is the minimun number of PMUs to be 
deployed in order to have the observability. 

A second issue is to find an “optimal” placement. Now 
we propose a simple, even if onerous, strategy for placement. 
After hnding a minimal number of sensors and a place that 
guarantee the satisfaction of Lemma we propose a greedy 
PMU placement procedure based on the sequential addition 
of one PMU at a time able to provide the best performance 
improvement, verihed by Monte Carlo simulations. For every 
possible new place for PMU, we run num run Monte Carlo 
simulation, of length TSTOP. For each of them we choose, 
randomly, the initial a, the switch i that changes its status 
and the time r of the action. The place for the new PMU that 
performs the minimum number of errors is then chosen. 

VI. Switching action in the non-ideal scenario 

So far, we considered the case in which the measurements 
devices were not affected by noise and the loads were static, 
which is not a realistic case. 

The measurement apparatus in a bus where measurements 
are taken is formed by a PMU and by a potential transformer 
(PT). A PT for metering aims to reduce the voltage magnitude, 
in order to make it measurable by a PMU. Both the PMU and 
the PT introduce errors to the voltage measured. We model the 
output of our PMU placed at bus j at time t by 

Vi (f) = (f) + G (^) + bj (t) (33) 

where ej S C is the error introduced by the PMU, while 
bj G C is the errors introduced by the PT. A common 


Algorithm 2 Greedy PMU placement 

Require: V such that the network is observable, num run, 
TSTOP 
1: min ^ 00 
2: minimizer 0 

3: for every possible place j ^ V do 
4: eiT ^ 0 

5: for t = 1 to num run do 

6: choose a ~ Z^([0,1]’’) 

7: choose f ~ Z//([0,... ,r]) 

8: choose r ~ ..., TSTOP]) 

9: run Monte Carlo simulations 

10: err = number of errors 

II: if err < min then 

12: min ^ err 

13: minimizer ^ j 

14: end if 

15: end for 

16 : end for 

17: V ^ V U {minimizer} 


characterization of the error is the total vector error (TVE). 
For example if x is the variable to be measured, and is the 
measured value, the TVE is 

TVE = ^ (34) 

Eurthermore, the loads are not static but they have a natural 
dynamic. In this paper we assume that the loads have constant 
power factor, and consequently 

= 7j. Vf, j = 2,..., n. (35) 

V{t)j 

We model the active power and the reactive power consump¬ 
tion at each load by 


p{t + 1)^1 =p{t) -I +np{f) (36) 

q(t)-i = diag( 72 ,... , 7 „)p(f)_i (37) 

where np{f) is a Gaussian random vector. This could seems 


too rough and unsophisticated, but we will see in section VII 
that is actually enough accurate. 

If we take into account ([3^, ( |3^ and ( |J7| ), the trend vector 
becomes 


^(^ 1 ,^ 2 ) = Ip{X„(^ti) - Xa{t 2 )) 

+ ('ll ~ ^*2 + 


S{t2 


+ S-ti — 6(2-1- 


IvXcr(ti) 


Un 


Un 

ti-i 

^ np{t) - iuqit), (38) 

t=t2 


being = p(f 2 )-i + Since a PT is a 

passive device and since in the following we will consider 
only fast measurement rate, we approximate satisfactorily the 
errors introduced by a PT with a constant, i.e. we will assume 
that bj{t) = bj. As a consequence, the term bt^ — bt 2 in ( [38] l 
vanishes. Because of measurements noise and loads dynamic. 
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the trend vector is typically non-zero even if there has not 
been any switching action, and the projection (|2^ may be 
almost one, leading to false topology detection. When any 
switch is closed or opened, since we are adding or deleting 
a branch, we are changing the curi'ents flows, reflecting on an 
abrupt, greater voltages variation. Therefore, a first strategy 
to avoid false positive is to not consider trend vector whose 
norm is lower than a defined threshold, called in the following 
min_norm. Moreover the additive noise can make the max¬ 
imum projection value of the trend vector maxC onto the 
library vectors considerably lower than one, even if a topology 
change occutTed. This fact prompts us to use a value for 
min_proj lower than the one considered in Algorithm [T] Of 
course, the use of a lower threshold makes the algorithm more 
vulnerable to false positive. The following example will give 
the idea for a possible solution. 

Example 3: Assume the grid without load variation and 
measurements noise, and that at time ti the f-th switch change 
its status. Consider the trend vector 

t) = y{t) - y{t - t). 

For t < ti and t > ti + t the projections of the trend vector 
onto the library are all equal to zero, because 

S{t,t-T) = y{t) - - r) = 0 

Instead for ti <t <ti -\- t, the trend vector is 

t-T)= A"^^ 

leading to a cluster of algorithm time instant of length t (or 
j seconds), in which the maximum projection coefficient will 
be almost one. 

A possible solution is thus to consider a trend vector built 
using not two consecutive measures, but considering measures 
separated by r algorithm time istants , i.e. 

t) = y{t) - y{t - t). 

and to assume that a topology change has happened at time 
t when we have a cluster whose length length cluster is r 
of consecutive values of maxC greater than min_proj. This 
idea will be clarified with some simulations in Section IVllII 
The former observations lead to the Algorithm for topology 
detection with measurements noise and load variation. 


Algorithm 3 Topology transition detection with noise 

Require: At each time t, we are given the variables a{t — 1), 
minimizer(f — 1), length_cluster(f — 1) 

1: a{t) ^ a{t — 1) 

2: each PMU at each node j record voltage phasor measure¬ 
ments yj{t) 

3: the algorithm builds the trend vector 6{t,t — r) 

4: if \\6{t, f — t)|| < min_norm then 
5: 5{t, f — r) 0 

6: minimizer(f) = 0 

7: length_cluster(f) = 0 

8: else 

9: the algorithm projects 6{t,t — r) in the particular 

library obtaining the set of values 





10 : if maxC >min_proj then 

11: minimizer(t) = argminC 

12: if minimizer(f) = minimizer(f — 1) then 

13: length_cluster(<) -s— length_cluster(f — 1) H- 1 

14: if length_cluster(f) = t then 

15: a{t) ^ minimizer(f) 

16 : end if 

17: else 

18: length_cluster(f) ^ 1 

19: end if 

20 : end if 

21 : end if 
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VII. Dynamic characterization of actual loads 

Load behavior at the individual customer level and at 
high time resolution can be a critical question in distribution 
networks. Due to lack of accurate and high resolution measure¬ 
ments at meters, there is no a clear answer to that question. The 
most commonly available data from loads come from meters 
with hourly (or 15 minutes) time intervals. For applications that 
are based on network parameters in shorter time than hourly, 
large uncertainty is caused by low resolution load data. For 
our topology detection algorithm, the load variation is very 
important, since it affects the topology detection capability, as 
( |3^ shows. To add more practicality to the proposed topology 
detection algorithm, a load measurement data set for five 


Figure 1. Load variability and load change variability based on normalized 
load values. 


houses in the United States is used. Load demand (kW) is 
recorded every second for a week. Statistical analysis of these 
load data is presented in Table |I] Table [TI] reports instead the 
statistical analysis of the variation of the load demand between 
two consecutive seconds. The field “Relative SD” expresses 
the standard deviation of the load or of the load change as a 
percentage of the load peak value, written in Table Figure 
[T] shows the histogram of load duration and load changes. 
The orange bars represents measurements data points over the 
percentage of pick load values. The blue bars show load change 
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Table I. LOAD VALUES 



SD (kW) 

Max (kW) 

Relative SD (%) 

House 1 

0.415 

5.990 

6.93 

House 2 

1.309 

10.956 

11.95 

House 3 

1.955 

11.578 

16.88 

House 4 

2.566 

12.364 

20.75 

House 5 

1.309 

8.155 

16.06 

Aggregate 

4.204 

27.229 

15.44 


Table II. 

Load variation 


Mean (kW) 

SD (kW) 

Relative SD (%) 

House 1 

0.000 

0.045 

0.11 

House 2 

0.000 

0.070 

0.64 

House 3 

0.000 

0.113 

0.98 

House 4 

0.000 

0.110 

0.89 

House 5 

0.000 

0.046 

0.56 

Aggregate 

0.000 

0.184 

0.68 



in one second intervals over the percentage of pick load values. 
In case of load changes, more than 90 percent of measurement 
time points have loads with less that one percents of pick load. 
It means that during one second time intervals, there is not 
huge difference in load data. In the United States, a num¬ 
ber of houses are connected to one distribution transformer. 
Therefore, the aggregated hve houses loads are considered 
as the reference for load variability in this paper. The lower 
the measurements frequency, the higher the load variability. 
Resampling the data, we obtain the aggregated characterization 
reported in Table III for different measurements frequencies. 


i.e. for one measure every second, one measure every hve 
seconds and one measure every ten seconds. Further numerical 
analysis shows that the load differences, for all the frequencies, 
can be modeled as Gaussian random variables thus validating 
equations ( [36l ) and (|3^, that can be assumed as a realistic 
model of load variation. 


VIII. Results, Discussions and Conclusions 
We tested our algorithm for topology detection on the IEEE 
33-bus distribution test feeder ini, which is illustrated in the 
Eigure 1^ In this testbed, there are hve breakers (namely Si, 
S 2 , S 3 , S/i, S's) that can be opened or closed, thus leading 
to the set of 32 possible topologies Ti,..., T 32 . Because of 
the ratio between the number of buses and the number of 
switches, some very similar topologies can occur (for example 
the topology where only Si is closed and the one in which only 
S 2 is closed). In the IEEE33-bus test case, Assumption[2about 
line impedances does not hold, making the test condition more 
realistic. Each bus of the network represent an aggregate of hve 
houses, whose power demand is described by the statistical 
Gaussian model derived in Section |VII[ dependent on the 
sampling frequency, and whose characterization is reported 
in Table III Regarding the measurement noise, we assume 
that the buses are endowed with high precision devices, the 
pPMU ifTSl . and with PTs. PMUs measurements are effected 


Table III. LOAD VARIATION FOR DIFFERENTS FREQUENCY 

Mean (kW) SD (kW) Relative SD (%) 

/ = 1 Hz 0.000 0.184 0.68 

/ = 0.2 Hz 0.000 0.425 1.56 

/ = 0.1Hz 0.000 0.604 2.22 


Figure 2. Schematic representation of the IEEE33 buses distribution test case 
with the five switches 



I [sec] 


Figure 3. Trend vector norm with noiseless measurements and loads non 
time varying 


by Gaussian noise such that TVE < 0.05%, based on the 
PMU manufacturer test information. It is also comply the IEEE 
standard C37.118.1-2011 for PMUs CSl. The PTs instead 
introduce a constant error (but peculiar of each node) that 
satisfy the requirements of the standard EOl . 

A. Trend vectors and noise treatment strategy 

Here we provide simulations that show the problems given 
by noise and that validate the strategy we use to overcome 
these problems. The overall time window simulated is of 
1000 seconds, the measures frequency is of one every 10 
seconds and the topology transition (/ = 1/10 sec“^), from 
(Ti = (1,1,1,0,1) to (72 = (1,1,0, 0,1), happens at f = 480 
sec. In Eigurej^we see what happens to the trend vector norm, 
while in Eigure|^we plot maxC, when we are in a noiseless 
scenario with loads not time varying. We can see that the 
switch action instant is clear. In the case with measurements 
noise and time varying loads, the trajectory of the trend vector 
norm, compared with min_norm numerically setted to 0.05, 
is reported in Eigure In Eigure instead we plot maxC, 
without putting to zero the trend vectors whose norm is lower 
than min_norm before the projection. We see that there are 
several cluster of time istants in which maxC is greater than 
min_norm, and this leads to a number of faulse positive. In 
Eigure]^ we plot maxC, after putting to zero the trend vectors 
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Figure 4. Trajectory of the maxC with noiseless measurements and loads 
non time varying 


Figure 7. Trajectory of the maxC with noisy measurements and time varying 
loads after putting to zero the vectors whose norm is lower than min_norm 
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Figure 5. Trend vector norm with noisy measurements and time varying 
loads 


whose norm is lower than min_norm before the projection. We 
still have more than one cluster, but only one has a length of 
Tj seconds, thus revealing the switch action instant, and the 
validity of our strategy. 

B. Simulation of the algorithm 

Here we tested the entire switches monitoring algorithm, in 
different situations and with different placements of PMUs: 

• 'Pas = {1) • ■ •) 33}, where every node is endowed with 
a PMU; 

• Vis = {3,8,9,10,12,15,16,17,18,19,21,24,25,27,30} , 
where there are 15 PMUs, almost one every two nodes; 

• Vj — {9,12,15,18, 24, 27, 30}, where there are 7 
PMUs, almost one every four nodes 

Vi 5 and Vj has been computed using Algorithm The algo¬ 
rithm has been tested in each condition via 10000 Monte Carlo 
simulations.. Firstly we tested the algorithm in the ideal case 
without load variation nor measurement noise, resulting in no 
etTors. Therefore, we can see that in the steady-state condition 



Figure 6. Trajectory of the maxC with noisy measurements and time varying 
loads 


Table IV. SIMULATION RESULTS WITH 33 PMUs 


SD [kV] 


non 

wrong 

decision 

total 

perc. of 



detections 

detection 

errors 

errors 

errors (%) 

0 


0 

50 

50 

100 

1.00 

0.68, (/ = 

1 Hz) 

0 

64 

67 

131 

1.31 

1.56, (/ = 

0.2 Hz) 

17 

131 

152 

300 

3.00 

2.22, (/ = 

0.1 Hz) 

72 

211 

244 

527 

5.27 


and in the absence of noise, the algorithm is extremely efficient 
for the 33-bus test case. It also overcomes the linearization 
from Proposition [T] and the initial Assumption [T] Secondly 
we added the measurement noise and different levels of load 
variation (or alternatively of measures frequency). The results 
are reported in Table IV Table [V] and Table VI The field “non 
detection” refers to the number of run in which the algorithm 
doesn’t comprehend that there has been a switching action, 
the field “wrong detection” refers to the number of times the 
algorithm detects a false action, while “decision errors” is the 
number of times the algorithm provide a wrong breakers status 
estimation. Of course every time there is a wrong detection, 
we have a decision error, too. If we subtract the values of the 
second column to the values of the third, we find very small 
number, meaning that, once detected the exact action time, the 
algorithm works very well. The main challenge is therefore 
to detect the topology transition time. However our approach 
is very robust with the number of sensors; the percentage of 
errors with only 7 PMUs (where we have partial information 
on the grid state) is always worst than the one with 33 PMUs 
(where we have global information on the grid state) less than 
2.5%, and probably with the best parameter tuning imin_norm 
and min_proJ), it can still be improved. 


Table V. SIMULATION RESULTS WITH 15 PMUs 


Relative 


non 

wrong 

decision 

total 

perc. of 

SD (%) 


detections 

detection 

errors 

errors 

en'ors (%) 

0 


0 

52 

52 

104 

1.04 

0.68, (/ 

^ 1 Hz) 

0 

73 

76 

149 

1.49 

1.56, (/ 

= 0.2 Hz) 

29 

135 

158 

322 

3.22 

2.22, (/ 

^0.1 Hz) 

74 

213 

252 

539 

5.39 


Table VI. 

Simulation Results with 7 PMUs 


Relative 


non 

wrong 

decision 

total 

perc. of 

SD (%) 


detections 

detection 

errors 

errors 

errors (%) 

0 


0 

56 

56 

112 

1.12 

0.68, (/ 

^ 1 Hz) 

2 

180 

185 

367 

3.65 

1.56, (/ 

:= 0.2 Hz) 

31 

199 

209 

441 

4.41 

2.22, (/ 

= 0.1 Hz) 

76 

245 

298 

619 

6.19 
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IX. Conclusions 

In this paper we propose a novel strategy for the monitoring 
and the identification of switches action in a distribution grid. 
The novelty of this algorithm is the possibility of running it in 
real time and the fact that it could work satisfactorily with only 
a partial knowledge of the grid state. It allows us, exploiting a 
voltage phasorial measurements, to understand both the time 
of the switching action and the new topology, by comparing 
the trend vector, built by data, with other vectors contained 
in a library, that represent the a priori knowledge of the 
electrical network. The algorithm has been tested in a realistic 
scenario, where both measurement noise and load variation had 
realistic characterization. In particular, we provide a simply but 
plausible model of the fast time scale load variation, validated 
using real field load data. The simulations show the algorithm 
behavior in different scenario. In the ideal one, with static 
loads and perfect measurements devices, the algorithm gives no 
errors. When instead we have noisy PMUs and load variation, 
we still have satisfactory results, strengthened by the fact that 
the algorithm performance are very similar both in the case of 
one PMU per bus or of one PMU every four nodes. 
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